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VARlABLfi RATE CONTlNtJOUS MODE SATELLITE MODEM 

FIELD OF XEIE INVENTION 
The present iavention relates generally to conmiimicalions and more particularly 
relates to a variable rate continuous mode modem suitable for satellite communications. 

5 BACKGROUND OF THE INVENTION 

Recently there Im been a vvidespread ]>rolifetatic(n of! satebite ^mmuni<^tions 
systems tiiroughout ttie world. A large number of consumets have acquired satellite 
communications equipment for home rccq)tion of satellite broadcast television signals in 
addition to many businesses having acquired such equipment for Hoc purpose of receiving 
1 0 and/or transmitting useful business information. 

Typically, satellite communications systems have been constructed using a 
combination of RF oiwrnitry, IF and baseband circuitry, A common trend in tfie art is to 
implement the IF and baseband modules using one or mpre ^PGi^^s an^^or ASICs. A 
disadvantage of this is that it is difScult or even impossible to periform some tcohniques tfiat 
15 ate only possible to perfonn in software. Attempting to perfopn these techniques in circuitr^^ 
if possible, would result in large, complex and expensive silicisn. 

Thus it is desirable to have a satellite modem Therein the recover portion is mainly 
implemented in software so as to reduce cost, making necessary changes much simpler to 
perfomi and potentially reducing ihe size of the receiver. 
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SUMMARY OF THE INVENTION 
The present inveation overcomes the disadvwtages of the prior art by providing a 
satellite modem ^ose baseband modules can be impiemectted mostly in software runmng on 
a suitable processor such as a high $peed digital signal processor (DSP) which are commonly 
available today. The modem includes a RF portico and a baseband portion. The digitized 
output of a matched filter is input to the baseband processing portion that consists of software 
executing on a DSP. 

The modem comprises an antenna coupled to an upconverter/downconverter. The 
upconvertcr/downconverter is coupled to a KF tcanscdver which ftinctions to receive a 
transmit signal &om a baseband transmitter module and generate a signal for input to the 
upconvert^/downconverter, The RF transceiver also functiotts to receive a downoonverted 
signal and output an IF signal for input to an IF module. The IF ftmctions to generate I and Q 
outputs that are input to a digital baseband module. The baseband module can be 
implemearted as software executing on a DSP, 

A xiovel aspect of the present invention is the method of performing phase acquisition* 
hi particular, tlie invention disclosed a tedmique of modulation wipe off v&erein fbe efifects 
of modulation are removed. The technique involves rotating all fhc symbol vectors into a 
single quadrant This is done in order to avoid the case where the constellation points are 
spread over more than one quadrant due to excessive noise being present in the chazmeh 

Another novel aspect of the invention is the signal detection technique. The signal 
detection technique performs signal detection and frequency acquisition in the presence of 
multiple signals in accordance with a method described hereinbelow. 

Still another novel aspect of the present invention is the timing acquisition technique 
that includes an mwrapping technique. This functions to correct timing errors before 
averaging is performed. In addition^ linear changes to the time estimates axe performed using 
LSR techniques. 

The ^tellite modem of the present invention, wherein the baseband receiver and 
transmitter are mainly implemented in software, provides fiwr a reduction in cost, makes 
performing necessary changes much simpler and reduces the size of the modem. In addition, 
the invention represents a significant step toward the goal of a coniplete software radio. 

There is thus provided in accordance with the present invention a satellite modem, 
comprising a housing, a antenna for recei^^g and transmitting radio frequency (RF) signals, 
an upcottverter/downconverter coupled to the antenna and adapted to i?)convttt RF signals 
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ftom a first firquejicy band to a second frcqtieaoy band, and adqrted to downcouvert RF 
signals from tiie second ftequeacy band to the first ficqtumcy band, an RF teansceiver 
coupled to ib» Tqjcoiivffller/downconvcrter and adapted to receive an ou^nt transmit signal 
and to generate an cvitput RF signal Acrefiom, the RF transceiver adt^d to receive an inpm 

5 KF signal fiom the upcoxrvertei/doMmcoiiverter and generate an iapxit receive signal 
there&om, an intcnnediate fteqiioncy (IF) module adi5>ted to receive the input receive signal 
and generate I and Q signals in response thereto, a baseband module adapted to receive the I 
and Q signals and to gencscate receive data in accordance therewifli and a baseband module 
adapted to generate the ou^t transmit signal in accordance vfiih a transmit data signal iapiA 

10 thereto. 

There is also provided in accordance vAUx die present invaition a receiver baseband 
apparatus, oomptiang hiput means adiqptod to receive an I and Q signal, an I matched filtM- 
adapted to receive the I agnal and generate an I filtered output herefrom, a Q matched filter 
adapted to reodve the Q signai and generate a Q filtered output ihere&om, a processor 

15 programmed to perform automafio gain control (AGC) and generate an AG control signal 
theieftom, perform timing detection and geawate an A/D dock control signal therefrom, 
perform phase detection and generate a voltage controlled oscillator (VCO) control s^al 
therefrom, a decoder aiapted to receive Ihe I output signal and the Q oulput signal fiom the 
processor and to generate a decoded output ttterofiwm, a deinterleaver adapted to generate a 

20 deintedeaved output in accordance with the decoded output signal ii^ut thereto, a forward 
error correction decoder adapted to generate ou^ut receive data in accordance with the 
ddntCTleaved output signal input thereto and a controller adapted to manage and control the 
mput means, I matched filter, Q matohed filter, the processor, the decoder, the dMdtKleaver 
and the forvraid error correction decoder. 
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BRIEF DESCRIPTrON OF THE DRAWINGS 
The mvMtiott is herein described, by way of example only, witii reference to the 
aocompanying drawings^ wherein: 

Fig- 1 is a block diagram illustrating the RF, IF and baseband portions of the variable 
5 rate satellite modem of the pxesent invention; 

Fig. 2 is a block diagram illxistradng the baseband module portion of the variable rate 
satellite modem in more detail; 

Fig. 3 is a flow diagram illustrating the xnain soflvwe processing loop of tite variable 
rate modem; 

10 Fig. 4 is a flow diagram illustrating the porocessing performed upon power up or after 

the communication link is broken; 

Fig. S is a flow diagram illustrating the hdHsl automatic gain control processing 
portion ofihe variable rate modem; 

Fig. 6 is a chart illustrating the results of correlation between ibjs spectral estimation 
IS and the spectral mask of ^ expected signal with both desired and adjacent signals having bit 
rates equal to 9600 bps; 

Fig. 7 is a chart illustrating the correlation result with the desired signal having a bit 
rate equal to 9600 bps and an interfering signal havmg a bit rate of 12S Kbps; 

Fig« 8 is a flow diagram illustrating the signal detect and frequency acquisition 
20 processhig portion performed upon power up; 

Fig. 9 is a flow diagram illustrating the timing, AGC, firequency and phase acquisition 
processing portion of the variable rate mod^ 

Fig. 10 is a flow diagram illustrating the fine frequency estimation processing portion 
of the variable rate modem; 
25 Fig. 11 is a flow diagram illustrating tiie coarse phase acquisition processing portion 

of the variable rate modem; 

Fig. 12 is a flow diagram illustrating the fine phase acquisition processing portion of 
the variable rate modem; 

Fig. 13 is a block diagram illustrating tracking portion of the variable rate modem; 
30 Fig. 14 is a flow diagram illustrating flie pro-track processii^ portion of the variable 

rate modem; 

Fig. 15 is a flow diagram illustiating the AGC traddng portion of the variable rate 
modem; 
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Fig, 16 is a block diagtam illustiating a model of the timing teddng portion of the 
variable rate modem; 

Fig. 17 is a block diagram illustratiiig the phase tradkmg portion of the variable rate 
modem; and 

Fig- 18 is a flow ^i^sjr^rtt illustrating the trackmg processing portion of 4e variable 
rate modem. 
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DETAILED DESCKIPTION OF THE INVENTION 
Notation Used Throughout 



The foUoMving notation is used Huoug^ut this document 



Term 


Definition 


AGC 


Antomatic Gain Contiol 


ASIC 


AppUcatioxx Specific Integrate Circuit 


BER 


Bit Error Rate 


DDS 


Direct Digital Synthesizer 


DEC 


Ciicolar Drift Fmctioii of EC 


DFT 


Discrete Frequency Ttansfoim 


DSP 


Digital Signal Processor 


ES 


Estimated Spectrum 


FB 


Frequeticy Estibnate 


FFT 


Fast Foiirrier Traaasform 


FIR 


Finite Impulse Response 


FPGA 


Field Programmahle Gate Array 


IF 


Intermediate Frequ^y 


IIR 


Infinite Impulse Response 


LIM 


Limit of Uncertainty Range 


LPF 


Low Pass Filter 


LSR 


Least Squares Regression 


NCO 


Numerically Controlled Oscillator 


NFFT 


Number of FFT points 


PI 


Peak Indexes 


PLL 


Phase Lock Loop 


QPSK 


Quadrature Phase Shift Keying 


RF 


Radio Frequency 


RMS 


Root Mean Square 


SNR 


Signal to Noise Rado 


TS 


Theoretic Spectrum Mask 


VCO 


Voltage Controlled Oscillator 



5 General Description 

The present invention is a variable rate continuous mode satellite modem C£^able of 
transmitting and receiving at a symbol rate in the ran^ of 9.6 Kto 256 K symbols/sec. Note 
that the bit rate depends on the coding rate. For example, at a code rate of i4, with QPSK, the 
bit rate is equal to the symbol rate, At a code rate of with QPSK, the bit rate is 50% more 
10 than the symbol rate. To aid in understanding the principles of the present invention, the 
description of the processing performed in the modem is divided into five main portions: (1) 
the main loop including a state machine for defining acquisition, tracking, etc., (2) signal 
detection and initial ftequency estimation, (3) acquisition of timing, frequency and phase, (4) 
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pre-tracldxig (AGC, timing and phase trackbig) vdHi a phase loop lutvmg a larger bandwidth 
and phase ambiguily to be solved and (5) tracking of the AOC, timing and phase. Additional 
processing includes the soft dedsion data detection. 

The main loop will now be described in more detail. A block diagram illnstrating the 

5 RFj IF and baseband portions of the variable rate satellite modem of the present invention is 
shown in Figure L The ZBodwx, generally referenced 10, comprises a housing (not shown), 
an iq)converter/downconv^ter 13 coupled to m e>ctemal antenna 12 and a KF transceiver 14 
coupled to the upconverter/downconvcrter 13, The upconverter functions to convert the 
output of the KF transceiver 14, e.g., L band, to the RF output frequency, e»g,, C or Ku band. 

10 The downconverter performs downconversion of the signal received over the antenna from 
RF frequencies, e.g«, C or Ku band, to L band for input to the RF transceiver 14. The RF 
transceiver downconverts the received RF signal to an IF signal. The IF signal output of the 
RF transceiver 14 is input to the IF module 16 vAdch fhnotions to gen^te the I and Q 
baseband signals. The baseband module IS demodulates the I and Q signals to yield tibe 

IS receive Rx data. The baseband module IS performs signal processing tas^ on the I and Q 
data from the IP module 16 such as Vitexbi decoding, de-mtcrleaving and Reed Solomon 
decoding 

In the transmit direction, the baseband transmitter 20 generates a baseband transmit 
signal in response to the input transmit Tx data. The baseband transmit signal is input to the 
20 RF transceiver 14 which i^converts tixe signal to an RF frequency signal which is then input 
to the external ant^ma 12. 

Baseband Module 

A block diagram xUxistrating the baseband module portion of the variable rate satellite 
modem inmore detail is shown in Figure 2, The baseband module 18 receives t}ie baseband I 

25 and Q signals. The analog I signal is converted to digital by A/D converter 2L The 8-bit 
digitized I signal is input to a matched filter 23. The analog Q signal is converted to digital 
by A/D converter 22. The 8-bit digitiTed Q signal is input to a matched filter 24. The 
matched filths 23, 24 are implemented in hardware but may altemadvely be inxplem^ted in 
software. The I emd Q outputs of the rmatched filters 23, 24 are input to processor 25, 

30 In one embodiment of the invention, for «>cample, the RF receive frequency fer C 

band operation is approximately 4 GHz and the RF transmit frequency is B;pptoxknsXeiy 6 
GHz. The RF transceiver 14 dow)converfcs/upconverts to/from an IF frequency of 
approximately 900 MHz. The IF module 16 downconverts the IF signal to analog I and Q 
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baseband signals within Uie bandwidth range of ^approximately 0 to 170 KHz. The baseband 
module IS demodulates Ifae baseband I md Q signal to a soft decision digital data signal. 
Note that the IF frequency, between the RF transceiver 14 and the IF module 16, may be 70 
MHz ±1B MHz or 140 MHz ±36 MHz, The KF transceiver 14 may comprise a 
5 commeroially available transceiver suitable for operation with the satellite modem of the 
present invention. For Ru band operation, the RF receive frequency is approximately 12 
GHz and the RF transmit frequency is approximately 14 GHz, The RF transceiver 
downconverts/iqpconverts to/from an IF frequency within the range of acpproximately 950 to 
1450 MHz. 

10 In an altemative embodiment, the signal from the RF transceiver constitutes a first IF 

which is it^ut to a modified IF module that outpxxts a second IF signal. This second IF signal 
is converted to digital via an AJD converter and input to a digital converter that frmctions to 
generate I and Q signals. 

Tbe processor 25 perfozms the mq'oiity of the signal processing functions in the 

15 modem 10 such as activity detection, timing acquisition, frequency acquisition* phase 
acquisition and tracking. The processor 25 junctions to output 3-bit I and Q signals that are 
input to a Viterbi decoder 26. The outpirt of the Viterbl decoder 26 is input to a deinterleaver 
27, The output of the demterleaver 27 is input to a Reed Solomon decoder 28 which 
functions to output the receive Rx data, 

20 A controller 29 controls and manages the operation of the baseband module and clock 

circuitry 31 fimctions to provide the various clocks neected by the elements of &e baseband 
module IS. 

Main Processing Loop 

A flow diagram illustrating the main software processing loop of the variable rate 
25 modem is shown in Figure 3. Upon power (st^ 32) a scries of start up tasks are 
performed. These include activity detection (step 34), acquisition (step 36) and pre-tracking 
(step 37). The activity detection (step 34) performs signal detection, AGC and coarse 
frequency acquisition. In this step, the modCTi is listraing for a signal from the transmitter 
The frequency/tuning offset is reduced to a suitable level such that when the PLL is started, 
30 the errors are relatively small and the PLL is enable of looking onto the signal. During the 
activity detect phase, no mjatchcd filtering is performed. In addition, the freqxiency is variable 
with an approximately frequency ofl&et of 30 KHz at a low rate and with a bandwidlh of 
approximately 7 to 8 KHz. 



During acquisition (step 36), AGC is p«rft«mcd as well as timtog, fine frequency and 
phase acquisition- Next, during tfie pre-tracking phase (step 37), phase tracking acquisition is 
perfoimed wherein a phase estimate is perfoimed every four symbols. 

Once the start up phase is completed, tracking is performed (step 38). This step 
5 includes AGC, timing and phase tracking. The ahove steps are repeated when tracking is lost 
on the sigcial or after each time Ihe mod^ is powered i^. 

Activity Detection 

Each of the above steps is described in more detail hereinbelow beginning wilh the 

activity detection phase. A flow diagram illustrating the activity detect method (step 34 
10 Figure 3) is shown in Figure 4. The foUowing steps are perfiarnied upon power up or aft 

catastrophic event whereby the communication link is completely broken* The initial 

automatio gain control (AGC) step includes measuring the energy at the output of the A/D 

converters 21^22 (Figure 2), i,e„ at the ii^ut buffer (step 40). 

the decimation of the input signal is performed (step 42), Decunation by a 
15 ^ctor of 1> 2, 4 or 8 is performed, Signal detection and initial coarse frequency acquisition 

and estimation is perfomied (step 44). If signal activity is not detected^ the method loops 

back and repeats steps 40, 42, 44 until a signal is detected (step 46). 

Each of the steps of the activity detection method will now be desmbed in more 

detail beginniiig witii the initial AOC st^ 40. A flow diagram illustrating tiie initial 
20 automatic gain control (AGC) processing portion of the variable rate modem is shown in 

Figure S. First, N samples of the hiput signal are collected whereby K may be 64, 128 or 256 

(step 50), The unaltered samples are e^q^ressed as 

^n^h^JQn (1) 

The mean square energy Ems is then calculated using the following equation (step 51), 

25 (2) 

wh^ Go is the gain of the AGC and IGo is tixe inverse gain. The inverse gain IGa is set equal 
to 256 * Eflis for scaling the RMS of the A/D converter input at 1/16 maximum input voltage 
level (step 52). The value of Go is delmnined from a lookup table (step 54). The G© value is 
then written to tiie AGC hardware (step 56). Signal detection is then performed (step 58). 
30 Signal decimation (step 42 Figure 4) is performed before the actual processing in 

order to improve the frequency resolution at lower data rates. The resolution achieved is a 
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particuter firaction of Uie symbol rate. It is preferable to have between 4 and 10 samples per 
symbol, excqpt for ±0 case of very low data rates. This means that the signal is relatively 
vddeband iu the nonnalized frequency domain thereby enabling a reasonable resolution 
i?sdthout reqidxing very lot^g FFTs. The resoliatioa can be expressed as 

binwidih fJN 



Symbol Rate 



Rs N 



(3) 



N 



>^ is the sampling rate 

i;^ is tiie symbol rate 

wyj is the number of samples per symbol 

JSTis the number of points in the FFT 
The paramieters of the decimation for activity detection, also referred to as gear shifting, is 
presented below in Tabic L 



Table 1 



Symbol Rate 
(KHz) 


Sample RiKte 
atA/D(KHz) 


MskBV/ 
(KHz) 


Dedmation 
Factor 


Number of 
samples per 
symbol 


Noaxialized 

Cutoff 

Freoucncv 


Relative 
Resohition(%) 


9.6-12 


768 


38.1 


8 


8-10 


0.0496 


3.1-3.9 


12-16 


864 


41 


8 


6.75-9 


0.0475 


2.63-3.5 


16-24 


928 


46^ 


8 


4.83-7.25 


0.04978 


1.88-2.8 


24-32 


512 


51.6 


4 


4-5.33 


0.1008 


1.5-2 


32-64 


512-1024 


51.6 


4 


4 


0.1008 


1.5 


64-128 


512-1024 


73.2 


2 


4 


0.1430 


1.5 


128-256 


512-1024 


116.4 


1 


4 


0.2273 


1.5 



The maximum bandwidtih is determined by the bandwidth of the higher rate plus a 
maximum ftequency ofiEset of 30 KHz. 

The relative required resolution is determined from the limitations of the phase 
acquisition algorithm. The phase acquisition algoritimi operate on blocks of eight symbols, 
trying four hypoGieses for the axis, each rotated by 22,5 degrees from one anotiber. 

Therefore^ over cigfht symbols, the maximum frequency ofi&et causes less than 10 to 
20% of 22.5 degrees, Le., less -ftian 2.5 to 5 degrees phase shift This implies a resolution of 
much less than one degree per symbol. This corresponds to a resolution of less than 0,3% of 
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the symbol period. This resolution, however, is very difficult to obtain, thus a fine frequency 
acquisition is petfonned prior to the phase acquisition. This yields resolutions between 1 .5 to 
3% of the symbol period 

The decimation is performed by filtering udlizing half band filters, *\^ch results in 
5 decimation by a factor of two. In addition, these filters are relatively inexprasive to 
implement since (1) approximately half tiae coefficirats are zero, (2) the filter order is 
relatively low due to a wide transition band between the pass and stop band and (3) only 
every second output need be computed due to the decimation by twd The design of half 
band filters is described in more detail in '"Multiiate Systems and Filter Banks," P.P. 
10 Vaidyanathan, Prentice Hall 1993; "A Tiiok for ihe Design of FIR Half Band Filters," P.P. 
Vaidyanathan and T. Nguyen, IEEE Transactions on Circuits and Systems, Vol. CAS-34, pp. 
378-389, April 1989, both of v^ch are incorporated herein by reference. 

In the embodiment presented herein, a filter having a length of 9 fbr ihe first stage is 
sufScient, while a filter length between 11 and 13 is sufficient for the second stage (or the 
15 single stage when decimation by 2 is performed). 

The coefficients for the filters arc presented below in Tables 2 through 4, 
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Table 2: Coefficients For Decimation By Factor Of 8 Using Three Filters 



-0,0336 


0 


FUter#l 

0.2835 0^000 0.2835 


0 


-0.0336 


-0.0422 


0 


Fiher#2 
0.2902 0.5000 0.2902 


0 


-0.0422 


0.0131 


0 


Filter H3 : Colirains 1 through 7 
-0.0236 0 0.0463 


0 


-0.0949 


0 


0.3144 


FUtM #3: Columns 8 tbroi^b. 14 
0.5000 0.3144 0 


-0.0949 


0 


0.0463 


0 


Filter #3; Columns IS Ifarongh 19 
'0.0236 0 0.0131 






Table 3 : Coefficients For Dedmadon By Factor Of 4 Using Two Filters 


-0.0426 


0 


Filter #1 

0.2904 0.5000 0.2904 


0 


-0.0426 


0.0146 


0 


Filter #2: Columns 1 4rough 7 
-0.0246 0 0.0472 


0 


-0.0956 


0 


0.3147 


Filter #2: ColnimDS 8 through 14 
0.5000 0,3147 0 


-0.0956 


0 


0.0472 


0 


Filter ^: Columns 15 ibrougb 19 
-0.0246 0 0.0146 






Table 4: CoeffidLeots ForDedmation By Factor Of 2 Using One Filter 


0.0169 


0 


Fater #3: Columns 1 through 7 
-0.0691 0 0.3041 


0.5000 


0.3041 


0 


-0.0691 


Filter #3 : Cohmms 8 through 1 1 
0 0.0169 
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Signal Detection and Frequency Acquisition 
This section describes the method used by the variable rate modem for signal 
detection and ftequency acquisition in the presence of multiple signals in the spectrum. 
According to Intelsat standards, two signals can be placed at firequencies that are a distance of 
5 0>7 * (symbol^ratei + symboLratez) as measured from center frequency to center frequency 

and vvdiercby the — ^ levels (fiffer by at most 7 dB, 

The metihods of the present invmtion utilize FFT periodograms for spectral estimation 
followed by cyclic correlation with the spectral mask of the ejcpeoted signal in order to detect 
the center of the signal. Figures 6 and 7 arc charts illustrating the results of correlation 

10 between the spectral estimation and the spectral mask of the expected signal. The result of the 
correlation for the Intelsat scenario, for example, is shown in Figures 6 and 7. In both 
figures, the esqjected signal has a power level lower by 7 dB than adjacent signals and has a 
dataratcof9600bps. 

With reference to Figures 6 and 7, multiple peaks can be detected in the correlation* 

15 In contrast to the sin^Je signal scenario, tiic maxiinal peak is not the sole candidate for the 
cttrtral frequency of the objected signal. Preferably, peaks in the correlation should be 
examined to detemiine the potential locations of the desired signal. 

Spectrum Estimation 

A flow digram iUuslrating tiie signal detect and frequency acquisition proccssmg 
portion (step 44 Figure 4) perforaied upon power up is shown In Figure 8. ITie first step 
performed is to calculate the spectrum of the received signal (step 60). The estimated 
spectrum (ES) is represented as an array or size 256, The following pseudo code describes 
tins step. 

Spect ffli ;pytiniaii on Pseudo Code 

L reset tiae ES array to zesro; 

2. for j=0 to numjreps - 1 (niim_reps « 32 Bar low— ^ levels and 16 for regular operation) 

get z , a vector of 256 received signal samples, 2, = + JQa i = 0,.255 
£S' = jF5+lFFT(f)|^ 

endfor 

ES 

3. calculate ES » ; 

ntm^reps 
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4. FFrahlfkES: 

m/) = P^^'~^^^^*'"^^*"^^. ^eie ES' is the FFT shifted version of 
^ ^ \ES(i + 128),/ - 0 -s- 127 

JSS, tooughout tiie remainder of this section, howcvet, JES" is denoted as ES; 

5. remove DC using: 

2) 



Signal Presence Detmnination 
Next, the presence of a signal is determined using the fbllowiiag pseudo code (stqj 

^^f^^l Fr ^sence Dctermingtion P??h4o Code 

1. smooth the Estimated Spectrum xisiiig an 11 coefficient low pass filter (LPF) 
presoited below in Table 5; 



^55(0 » S ESii-J)h(j% i « 0 ^ 255 
where = 0 for f < 0 and i > 255 ] 

Table 5: Coefficients For Low Pass Filter 
""X^ h;^ El ~ho 51 h^^ "hi[ 

0.010 0,0249 0.0668 0,1249 0.1756 0.1957 0.1756 0.1249 0.0668 0.0249 0.010 



62), 



2. calculate RES, the resolution of one bin in the FFT 

_ SampHngFrequenqy 
NFFT 

where NFFT is defined as the nmnber of FFT points; 

3. calculate chkjmg, the range of indmes to be checked in the FFT 

J m&xi3y<10^, 0J' symbol _rate ^ 
cHk_r^^ce.t^-^ — ^ 

4. calculate max chk and min_chk 



5. calcula:be act ratio 
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max cbk 
act ratio » — ; ,~ 
" mi3i_chk 

If the value of actjratio is greater than 1 .75, it is determined that a sdgnal is present (step 64). 
Otherwise it is deteixnined &at no signal is present 

5 Correlation Calculation 

The next step is to perfoim a correlation calculation (step 66). The following 

definitions apply: 

ES: estimated spectmm 
TS: theoretic spectrum mask 
10 EC: an airay of size 256 representing the coirelation fhnction between the non 

smoothed ES and the TS 

JJM: the limit of the range of mcertainty 

The correlation is detrained using the following pseudo code (step 66). 
15 CQi;yelatiQtt Detem;anfltiQp Pseudo Code 

1 . calculate LIM using the foUovdng equation; 

L JtES J 

2. calculate the following; 

20 £^5((jfc + f)^i,)rs(kXl2B -'Lim<i< 128 + Lim 

Peak Determination 

The next step is to perform a peak detemaination algorithm (step 68). The follo\ving 

definitions ^ply: 

25 DEC: circular difiEerenoe function of EC 

PI: list of peak indexes 
NJWIN: the minimal distance between peaks 



30 



The i^ak determination is determined using the following pseudo code. 



Peak Dete r^inaHnn Pseudo Code 



1. calculate DEC using tiie following equation; 

^ \EC{i) - MCii -l),tm 1.255 
35 2, find PI, for the list of iadexes / wMoh fulfill the conditions i e PI 



P-1318.US 15 



jDECQ) ^OwdDMCii^l) < 0,i = 0..254 
\DECi2SS) > 0 9iidDEC(o) < 0,f = 255 

3. find the valise of the minimum peak and eliinkaiate from tbe list those peaks that are 
10 dB (0.1) smaller tiban the maximiixn peak; 

4. calculate N_WIN using the following; 



N WIN ' 



floor(o35 ^^^ symbol rate ^64KHz 
\ nss J ^ 

Jloor^A75 symbol _rate>64KHz 



5. If two or more peaks are closer than NJWIN, keep the highest between them in the 
list. The order of the removal is from the highest peak to the lowest, Note that it is well 
known that the center frequencies of two signals arc separated by the sum of half the signal's 
bandwidth, i-e., QJisymbolj^ej-^symboljriXtei). In practice, smce (1) flae peaks of 
10 inteifewars are not accurately detected, (2) the peak of a small rate signal tends to get closer to 
the center of the high rate signal, and (3) the rate of the interferer is not known, iMs range 
should be decreased to the value given above. 



15 apply: 



Symbol Rate Calculation 
The next step is to calculate the symbol rate (step 70). The following definitions 

BW2: one half of the bandwidth of the signal 

The symbol rate is calculated utilizing the following pseudo code. 

20 Symhol M^x^^ Calculation Pseudo Code 

1 . calculate B W2 using the following; 

2. calculate the following; 

25 iVb_/eve? = JlfM7i({Ey(0|0.85-inm_c»* ^£5(0^1.^^ 

3. for /= 1 to the length of PI: 

3.1 define Area as: 

Area = ^\PI(l) -BW2<i< PZil) + BW2 

3.2 define Aiea2 as; 



30 Armi^U 



Pm - - BW2 < i < PI(l) +IbW2 

a o 



3.3 find max_val: 

max val= max ESJk) 
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3.4 redefine the indexes of Area2 as k such that ES^ (*) < 

Area2 = |Jp/(/) - 1 5r 2 < i < /-/(O + 1 5 ES, (0 < 

L I o 8 2 

3.5 calculate mean_val using; 

fnean_val = is^s^ 

ininiber of indexes in Area2 

3 .6 find Fn, the smallest index in Area greater than PiaX such Ihat 
^Ss(F^) ^0.5- (mean _ v<j/ + iVb _ /ev«/) 

3.7 find R, the highest index in Area smaller than PIG), such that 
^SsiF,) ^ 0.5 ♦ imean_val + Nb^level) 

u , if both Fu and FI are found thca"ialculate the estimated symbol rate 
(es_symboI_ratc) corresponding to ¥1(1) usingthc foflowing: somawa symooi rate 

es_symboi__rate = (JP„ -Fj)-ItES 
4. end for 

15 «c«ecJ;^? ^u*** P**^ «_-?);>»6o/_/-flre higher than 1.2 times the 

15 esqpected symbol rate or lower than 0.75 times the expetitsd sj^bol 

IftheUstofpeaks is exhausted, the method continues witii step 60 (step 72). Ifthe 
list of peaks is not exhausted, the peak closest to the center is selected (step 74). 

Frequency Acquisition 
Hxe next step is to perform firequency acquisition (step 76). The foUowing definitions 



apply: 



EE; the nonnal fiequen^ estimation 

pi: the peak index where a candidate signal is centered 



The ftequenoy acquisition is perfiarmed utilizmg the fi)Ilowing pseudo code. 
Frequency Acouisitinn Psaii^^ r^^^ 

1 . file niTOU fejpiency estimation (FE) is calculated using the following 

rr, = yjj — 128) • RES 

fi^queiy^ bS^^SS:^"^ and fi. central fiequ«.cy is shified in o^er that 

Conununlcations is then attempted (step 78). If communications are not successful, 
35 the particular peak is deleted from the list and the method continues wim Step 72). If 
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commumcation is successful, the meawd tctmlnates. Timmg acquisition and fine frequency 
acquisition are then perfoimcd (described in mote detail below), 



Ti ming, Frequency and Phase Acquisition 
A flow diagram illustrating the timing, AGC, frequency and phase acquisitioa 
5 processing portion (step 36 Figure 3) of the variable rale modem is shown in Figure 9. A 
matched filter begms operation at a nonnal gear shift at 4 samples per symbol (step 90). The 
matched filter is implemented in hardware and operate without a phase rotation. 
Alternatively, ti» matched filter can be implemented in software, The matched filter 
hardware controls die clocking of the A/D converters. 
10 The normal gear shift ratios are presented below in Table 6. 

Table 6: Normal Gear Shift Table 



Index 


LowSyaxboI 
Katt(KHz) 


Symbol Rnte 


Sainples/ 
Symbol 


Low 

Frccjuency 
(KHz) 


High 

Frequency 
(KHz) 


7 


9.6 


16 


64 


614.4 


1024 


6 


16 


32 


32 


512 


1024 


5 


32 


42 


24 


768 


1008 


4 


43 


64 


16 


672 


1024 


3 


64 


SS 


12 


768 


1020 


2 


85 


238 


8 


680 


1024 


. 1 


128 


170 


6 


768 


1020 


0 


170 


256 


4 


680 


1024 



First, 128 samples are coUeoted. At 4 seniles per symbol, tWs is equivalent to 32 
15 symbols. The first AGC acquisition is peifonned on ttie block of 128 sanies (step 92). 
Ehning this step, tiie fiillowing is calculated. 

" f-l 

where N - 128 samples at 4 samples per symbol and the AGC value is calculated 
utilizing the following equation: 

20 AQC^ ^'^ 



2x16* 

The next step is to perfijim timing acquisition (step 94). The foUowhig definitions 

apply: 

'block of symbols': 16 contiguous symbols 



til DFT estimate of a block of symbols, i « 1„N 
N: number of DFT estimates 
*group of DFT estimates': grox^ of N DFT estimates 
Ti: aveiage estimate of group of N DFT estimates 
5 K: number of average estimates Ti, which is the number of points used as input to the 

Least Squares Fit 

The timing acquisition is perfoimed iitilisdng the following pseudo code. 

In connection with the following Pseudo Code^ tfa^ foUowing assumptions are made. 
10 (1) blocks of 16 contiguous symbols at 4 samples per symbol are available, I.e., 16 * 4 = 64 
continxious samples; (2) the next block of 16 symbols arrive after a break of 16 * 3 = 48 
symbols; and (3) the maximum drift betweeix the receiver and the transmitter clocks is 
approximately lOOppm. 

The method of the invention utilizses the timing recovery technique described in M. 
15 Oerder and H-Meyr, 'THgital Filter and Square Timing Recovery," IEEE Transactions on 
Communications, Vol. COM-36, pages 605 to 612, May 1988, incorporated herein by 
reference. The method of the indention is given briefly below. 

1. Perform steps 2 through 5 for K groiqss made up of N DFT estimates, each 
calculated from blocks of 16 symbols. 

20 2. Calculate timing estimate % based on the DFT tedxoique for 16 contiguous 

symbols, i « 1 ,«,N, thereby obtaining N DFT estimates each based on 1 6 symbols. 

3. Make a histogram (3 bins in tiie range of -0»5, O.S) of the N DFT esthnates tj 
obtained in step 2. 

4. Classify the timing range (0, Va, Vi, the group the N estimates are in. The 
25 decision is (krived fiom logic based on the histogram in step 3. 

5. Unwrap N DFT estimates baaed on a range decision derived j&om logic and 
calculate Iheir average T;. 

6. Unwrap the K average estimates Ti and do a least square fit to the K averages to 
produce tiie final estimate. 

30 

Timing Acquisition Pseudo Code fin more detail^ 

L for k = 1 to K perform steps 2 through 5 inclusive )groups of average estimates T^); 

2. calculate the timing estimate ti based on a DFT algorithm for 16 contiguous 
35 symbols, i 1..N, thereby obtaining N DFT estimates each based on 16 symbols contiguous 

at 4 samples per symbol with no discontinuity in time. 

fori=l toN 
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'■r = S W + ^ • k W + 0 + 1 6 ■ 4 • p - ]- log^ + <S5SW • |z(4y + 2 + 16 • 4 • - 1))!'} 
5 {- log[l + 5aV« • jz(4i + 1 + 16 • 4 • (/ - l))f ]+ logjl + jSyW . (z(47 + 3 + 16 . 4 • 0 - 1^^^^ 



eiidfor(i= 1 toN) 
wliere 



^0 



and Es reprints energy per symbol. 
The prelimixiary estimate is fhcsa given by the following 

1 fe^ 
t^- — ^arctan — /«l»Ar DFT estimates 
2'pi \ej 

Note that the arctan fvEncdon can be approximated using Ihe least squaies method and 
the logarithm is the natural logarithm to the base *c\ 

3- construct a histogram of the N DFT estimates ti wherein 3 bins m the range -0.5 to 
0.5. 

3,1 count the nxmiber of estimates ^ in each bin as defined by 
bin 1 = [-0.5, -1/6] 
bin 2 = [-1/6, 1/6] 
15 bin 3 « [1/6, 0.5] 

32 sort counts of the three bins in descending ordet aod save their respective indices 
as follows: 

max = value of tiie maximmn eoxmt 
20 mid « value of the middle oonnt 

min = value of the minimum couot 

Imax - todex of the bin wife the maximum count 
Imid ^ index of the brnwilhi the middle count 
25 Imin index of the bin with the minimum coimt 

4. classify the portion of the timing range the group of N estimates are in; 

rauge 1 = around zero, Ic-, [-1/8, 1/8] 
30 range 2 « aroiind i.e.> [1/8, 3/8] 

range 3 = around 'A or i.e., [-'A, 3/8] and [3/8, Vz] 
range 4 " aroimd -i<s i-e,, [3/8, -1/8] 

if (max ^ 1 0) & (mid + min ^ 6) then 
35 ifW===2fhen 

range -1; time round 0 
elseif Imtot " 1 

range = 4; time negative 

else 
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else 



end 



raiige«2; time positive 

end 

if Imin ^ 2 & (max - mid ^ 4) 
iange'=3; time round S4 

elseiflttiia= 1 

OTagc = 2; lime positive 

else 

lange == 4; lime negative 



5. unwrap N DFT estimates based on tiie range decision axid calculate their average Tk 
using the following pseudo code; 

Average Cg^ culation Psefudo Code 

S > 1 . initialize thresholds: 

thimid=0.25 

tiirhlf-'O.OS 
5.2, calculate 



meam 



t , % av€3:age of DFT estimates before unwrapping 



5.3. do for group of estimates i = 1 to N 

sgn - 0 % case on value of the variable 'range' 
case 1 
NOP 
case2 

if ti < -thrmid then sgn = 1 ; 
c^e3 
if mean>0 

ifli<-thrhlfifaensgn=l; 
else 

if ti> liixWf then sgn« -1; 
esnd 
case 4 

if ti > thmiid then sgn = -1; 
end; 

end do for group of estimates 

5.4, calculate the average estimate of the group of N DFT estimates as follows 



endfor (k = I to K) 



Tj^ I = 1.JV",% calculate average cstunate 



6. unwrap the K average estttnates Ti and do a least square fit to the K averages to 
yield die final estimate; 

We assume that Ti depends linearly on time. The obfi^rvation intedrval, however, is 
limited to tibe range [-0.5, 0.5], 
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6. 1 unwrapping the K average estimates 

The first stop is to vawiap the K average estimates and receive new values Ti for 
i = for odd K. 

T^ywTftjipnp; P sgudo Code 

1. Detect positive and negative jumps in timing that are larger than 0.5. A positive 
jump is detected If 

TUp ^l) -TOP) > 0^5 
A negative jump is detected if 

rc;n+i)-rc/n)<o.5 

^where TQ) are the time estunatcs and j is the time index. 

2. For positive jumps at index jp, subtract 1 from all subsequent time estimates (with 
indices S: jp, 

3. For negative jumps at index jn, add 1 to all subsequent time estimates (wrtii indices 

2rjn, 

6.2 Calculate estimates of the slope a and the avetage b of the i 1 to usmg a 
least square fit estimation 

The estimate is calculated to be: 

2 

with no prediction, or 




with a prediction of p=l/2 orp«l. 
Where b and a are calculated by a least square fit to a linear curve given by 

e— (JC-l)/3 

and 

< W r 

for odd K, The indices are chosen so that flie first Ti generated is multiplied by 
-(K-l)/2. 
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With reference to Figure 9, once the timing has been acquired, the xuesA step is to 
apply a matched filter (step 96), H^te matched filter is applied using 1 sample per symbol at 
the filter oulput ^th no phase rotation. The next step is a second AGC acquisition on a block 
of 128 samples at 1 sample per symbol (step 98). 

At 1 sample per symbol, 32 samples are collected. The second AGC acquisition is 
pexfomied on the block of 32 samples. During this step, the followbcifi is calculated, 

ip^ere N « 128 samples at 4 samples per symbol and the AOC value is calculated 
utilizing the following equation: 

Ems 



10 AGC:^' 



2x16^ 



The next step is to perform fine frequency estimation (step 100). The fine frequency 
estimation is peirtormed on the output of the matched filter after timing acquisition has been 
perfonned. Note Ifaat in order to increase pcrfottnance at low SNRs, preferably two 
15 periodograms are averaged raftier than using only one. A flow diagram iUnstrating the fine 
fiequmcy estunadon processing portion of the variable rate modem is shown in Figure 1 0. 

The following definitions apply: 

NFFT: the number of bins in the EFT, 256 for high SNR and 1024 for low SNR 
PADD; zero padding fraction, 0 for high SNR and 0.5 for low SNR 



With reference to Figure 10, liie timing acquisition is performed utilizing the 
following pseudo code. 

Fine Preqnencv Acquisition Pseudo Code 



1, fori«OtoNFFT*PADD - 1 

collect the output of the matched filter (step 110) ^jQj at tljie correct 
time epoch (step 112); 

do Zi = zf in two stages: 
30 calculate zf ^u^+ jW; , wherei*^ = If - , - 2I^Q^ 

calculate = + /y,-, wherex. - w? ^'wf > J'i = 2tt,>v, 

2. calculate a NFFT point of {z^ with NFFT*PADD zeros padded (step 114), The 
absolute value of the FFT squared is denoted by ES. 

35 FFT shift ES: 
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NFFT _^ 



2 

where ES' is Ihe FFT shifted vision of ES. From here on in, £"5" is denoted by ES. 

3. find &c majdmum erf {eS}^'^ (step 1 16); 

4. £aid ^ index of tiie peak denoted by peal(;_index; 

5. the fiiequency estimation (step 118) is given as 



1 symbol rate ( r . j. NFFT^ 
4 A^^l^^^-^^-^J 



6. alternatively, better results can be achieved if the modulation wipe off is perfonned 
as follows 

7, m approximatioii of the modiiiatioii ivipe off can be achieved using the following: 

7.1 t>-Liog.(/f-^e?)J 

This can be perfonned using iterated comparisons. The iteration is through fbs values 

7.2 The wipe off is performed by 

2" 



With reference to Figure 9, the next step is to perform phase acquisition (step 102). 
The input to the phaw acquisition portion of fte method is Zn which a» complex samples Q 
andQ)^atatcofl sample per second, The value of N is equal to 16 complex samples. One 

assumption made is that the fiequency offeets are smaller than — - — where T is the 

symbol period, i.e., for N = 16, less than 1/128 of flic symbol rate that is approximately the 
maximum firequencyr error. 

The phase acquisition process receives N complrac samples (N = 16) and functions to 
detcimine the rotation angle ibsst brings tiie samples after modulalion -wiping to a minimum 
variance. I and Q are tiie in phase and quadrature components, respectively. 

^»**^-+/a.»-0..iV-l 

The mefliod consists of two stages: (1) a coarse grained search that chooses between 4 
•hypothesis' of signal rotation with modulation unwipkg in the first quadrant and (2) a fine 
estimation of tiie phase rotation angle. The role of the operation is to wipe off tije effect of 
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the QPSK modvdatioii on the phase. Before ithis operation, the phase of z is the sum of an 
unknown phase and of the phase due to modi^ation. Afl»t this op^adon, the phase dq>aids 
less on the data. 

A flow diagram illustratmg the coazse phase acquisition proccssmg portion of the 
5 variable rate modem is shown in Figure 11. A flow diagram illustrating the j5ne phase 
acquisition processing portion of the var^le ^te modem is shown in Figuic 12. 

With reference to Figures 11 aod 12, the pseudo code for the phase acquisition 
follows. 

10 Phase Acquisition Pseudo Code 

1 . Coarse Grained Search 

1.1 for hypothesis Hk, fc«-l, 2, 3, 4,...rotate Zn by adding the phase ejt k-22.5** (step 
120) in accordance wilh 

15 ^,(©)-2^-e^' 

1 2 wipe oSsnf^i^ modulation as £>llows (step 122); 
foreachz,«/^+/a 

check in whidi quadrant.it appears 
iflnOn < 0 then swap J„ and Qn 
20 take the absolute value pf and Qn such that is now in the first 

quadrant; 

endfor 

1.3 sum the wiped off (or rotated) vectors Zn(0i^ using the following (step 124) (note 
25 that the sum should continue until N-1 « 15) 

1.4 find the energy for each of the hypothesis (step 126) using 

At this point, it is relatively easy to sdbbow that for a rotational transfonnation^ the 
30 miinimurrji variance condition is equivalent to finding the maximum energy. 

1.5 choose the hypothes^ withmaximum energy (step 128). 

2. *Fine' Phase Estimate 

2.1 let S{&^)^X+ JY (step 130) be the vector having maxunum energy; if Y> X 
35 (step 132) then swap X and Y (step 134) and set a Jlag (step 136); 

2.2 calculate a = atan(Y/3C) ushig the techniques described below in 2.2.1 to 2,2.4 
(step 130); 

2.2.1 find a quadratic fit for the inverse'term u-l/x: 
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137); 



operations as foS^^' X to lie range of [0.5,.l] which can be perfonned u«ng left shift 

3c = 2* - a: 

K = 4.2349x' -5.83543c +2.6192 

22JZ calculate 

v-i^ = 2-.(y.«) 

the 'unnonnalization' by a factor of 2-« is performed using right shifts (step 

2.2.3 calculate a = atan(v) using a quadratic fit (step 138); 
a = -0.0052V* + 1.0804V -0.2886 

Note: (1) these coefficients are scaled in radians; 

(2) ^^^proximation assumes that 0<v'<l which is tme because v=x/y and 
2»i assures tb3.t X > y; 

2.2.4 if the swapjlaz « set (step 139) than a = 90° - a (step 140); 

given a^'^ '"^'^ ®' ""^^ subtracted fiom the phase of the samples is 

0 - ®j +a-45* (step 142) 
2.4 the rotation by ® yields 



This rotation operation brings in the received signal at phase ±fl:/4^«-/4, etc. in 
accordance wift fte data (for noiseless signals). Note that steps 132 through Z40 comprise 
the detail implementation of step ISO. 

The rotation angle ® obtained during acquisition is the initial phase of the phase 
control loop. 

A block diagram iUustrating pie-traddng and trackmg portion of the variable rate 
modem is shown in Figure 13. lids block diagram describes both the pre-traddng and 
tracking stages ofthe main loop processing. Hie I and Q data output ftom the matched filter 
isinputtoamixer(multipKer)211. The I and Q output of the mixer 211 is input to AGC 
estimation 213, timing detector 214 and phase detector 215, 

The AGC estimation 213 Junctions to generate an AGC control based on the energy 
or average power of the input signal. The AGC control is passed through a loop filter 216 
before behig ^lied to the AGC hardware hi the IF module section 1 6 (Figure 1). 

The timing detector 214 functions to generate a noisy timing estimate which is first 
filtered via loop filter 217 before being applied to the dock generation drcuitry (not shown) 
inihe A/D converter section ofthe baseband module 18. The A/D dock control signal output 
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firom the loop filter 217 functions to adjiist and correct for clock tixning ezrors in the A/D 
converters. 

The phase detector 215 functionf^ to g^erate a misy phase estimate that is input to 
the loop filter 218, One signal output of the loop filter comprises a VCO control signal which 

5 is input to the VCO circuitry (not shown) in the IF module 16. The VCO control signal 
fbnctions to rai$e or lower the frequency ou^ut of the VCO. A second signal output of the 
loop filter 218 is the phase estimate ^ 4at is fed into the e^^*^ block 212, This block 
functions to coxrect for phase and small frequency oSfsets. The output of the block 
constitutes the second input to the I/Q mixer 21 1 . 

10 The I and Q oxstpnt of the mixer 211 is also input to the soft decision block 219, The 

matched filter ou:^uts 5^, after time synchronizatiou> fi^quency and phase concction, must 
pass the soft decision block 219 before being input to Viterbi detection in the Viterbi decoder 
(not shown) in the baseband module 18 (Figure 1) (see also step 158 Figure 14). This process 
is performed on 1^ and Qn indcpcndentiy. The Viterbi decoder receives for each In and Qn, 3 

15 or 4 Wts of quantiatiott. The exact number depends on the specific Viterbi detector. One of 
the bits is the sign and the other bits describe the signal level. 

For es»cample, if the le^el thresholds are ai < a2 < as, the two bits of soft decision 
which describe the level are shown below in Table 7: 



Table? 


k 


Bit Repxcsentsdon 




00 


ai<lB<a2 


01 


ai < In < o.% 


10 


Ct%<lo 


u 



With reference to Figure 3, once the acquisition stage is complete, pre-tracking can be 
perfonned (step 37). The pre-tiack stage is in principle normal tracking of AOC and timing, 
while tiie phase tracking loop has a larger bandwidth and the 90° ambiguity still must be 
25 resolved. A flow diagram illustrating the pre-track processing portion of the variable rate 
modem is shown in. Figure 14. 

The method comprises (1) AGC tracking as described hcreinabovej (2) trming 
tracking using a method that opiates on a group of symbols (at one sample p^ symbol) to 
generate an estimate that is filtered using a smgle pole IIR filter; the method is called at a rate 
30 of approximately 160 Hz which 32 symbols being processed each time, and (3) phase 
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tracking performed once per symbol for symbol rates below 64 KHz aad once every 2 or 4 
symbols for higher symbol rates, llic criteria for declaring tracking loc^ 
high raw BER or one or mote lock detectors. 

Fixst^thoIaMQdataisiiipixtto&einatdiedfi^ AGC Hacking is then 

5 performed on the ou^ut of the matched filter (step 152). A flow diagram illnstratiDg the 
AOC tracking portion of the variable rate modem is shown in Figure 15. The pseudo code 
for the AOC tracking follows. 

AGC Tiap Kinff 

10 L initialize variables inchiding *coxmt' and set the threshold to 256 (stq) 180); 

2. check if connt= 8 (step 182); 

3. every T milliseconds (T «e 100 ms) receive the output Zn (n represents one sample 
per symbol) of the matched filter (step 186) while count is < 8; 

4. calculate ^„ as follows (step 196) 

15 g,^I^hold^\2,\' 

5. calculate AGC value (step 198). Filter g„ using a Buttcrworfh low pass filter (LPF) 
transfisr function shown below 

0- 073 'O+z'^) 

1- -0-854-Z- 



with a 3dB cutoff rate at 0.1 of the digital Butterworfh LPF. The corresponding 
20 <Uf&rence equation is as follows 

(?^,=0.854-i?,+0.073.(ff,+?^,) 

6. the ou^ut Yn of this filter is the new ^inverse AGC* (m lie energy domain) (step 

200). 

25 7. to prevent cHpping atthe input to the A/D converter, the threshold is adapted every 

K*Tms(Ki«8, 16). 

8. find i« in an i«v«r^e_^GC_toWe and write liie valiie ^ 
module; 

8. increment count by one; 
30 9. after N samples (N 8 in &e above example) xeset count to 0 and calculate Bxos 

(step 184); 

10. if Ems > 0.25 (step 188) then increase the thr^hold (step 192); 
Threshold^ a 'Threshold^ of « 1 + 1/8 

35 else if Ihe threshold > 1/256 AND Ems < 0.24 (step 190) then decrease threshold 

(step 194); 

Threshold - J5 - Threshold, « 1-1/32 

endif 



P-151g-US 



28 



With reierenoc to Figure 14, the next step is to perfenn timing traddng (step 154). 
The timing tracking portifln comprises a loop filter, an idfeal Jategrator and an inherent 
mimerically conttoUed osdllater (NCO). A block diagram illiatrating a physical model of 
the timing tracking portion of the variable rate modem is shown in Figure 16. The filtered I 
and Q data at 1 sample per symbol, i.e., the output of the mixer 21 1 (Figure 13), is input to a 
Mueller error detector 170. The oulput of Hie error detector 170 is input to a loop filter 172. 
The oulpot of the filter constitutes the A/D clock control signal, nus signal is uqmt to NCO 
174 that fbnctions to generate the clock for ihe A/D converter. The oirtput of the NCO is also 
fed back to the Mueller error detector 170. 

The timing signal modtel of the loop filter 172 can be represented as flie following 



and the NCO can be represented as 



Note that the loop filter has a gain of p, an ideal mtegrator and a zero at y, i.e., proportional 
and derivative contioL The timing tracking loop utilizes the Mueller and Mueller timing 
error detector 170 which operates on the filtered in-phase and quadrature components on a 
block of 32 symbols at a rate of qjproximately 320 Hz (deteraiined by tiie operating system). 
The timing error detector is given by tixe av«j«e of the I and Q error components: 

2K 

'Miiere 

=7, .sgnC/^,)-/,^ •sgn</,) 

and K * 1/16 is the nominal value of the signal amplitude as set by the AGC. Next, the 
average of a block size N 32 is computed as follows: 

The error process is filtered by the loop filter 172 as follows: 
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where is liie correctiott writtea to the hardware, Tlie NCO 174 or ideal totegrator is part of 
Ihe system and may or may not be implemented in softwiate. It should, however, be taken 
into account when calculating the dosed loop transfer fmctiotL Based on the timing signal 
model given above and the assomptioii that the error detector is ideal, i.ei3 it yields the exact 
S timizxg error), the closed loop transfer function of the estimated timing epodi is given by 

Note that it is convenimt to construct the filter using the classical parameters fiom 
continuoxis time in terms of ^2^^, the natural fiequency and 5, the dampitjg fector, it is 

desirable to express p and y in terms of o&n and |. This can be achieved in ttvo stages. The 
1 0 jStst stage is to convert to the discrete time domam standard form by solving Ihe following s et 
of equations: 

&om which wc obtain 

15 or for the other direction 

'y8 = 2(l-rcos^) 

"2(l-rcose) 

We now express r and S in terms of ^ design parameters <0n and § as follows. 

20 This is described in more detail ha Best, Phased Locked Loops, Second Edition, p, 342, 
McGaw Hill, 1993. In addition, a more detailed description of the timing error detector 170 
is given in EBEE Transactions on Communicadons, May 1976, pp. 516 to 530. 

The following is based on the assumption that iht timing error detector gain is unity. 
The defector gain will now be discussed in more detail. The transformation is based on the 

25 impulse invariance conversion ftom continuoiis to discrete time. From r and we can 
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compute and y. Based on the above equations, assuming an operatiog rate of 1/T - 320 or 
640 operations per second and setting ^ ~ 0.707^ the pazameters for p and y are presented 
below. 



Table 8.' Parameters >9 and y 


^(Hz)@320Hz 


/? 


r 


/„(Hz)@640Hz 


0.1 


0.0028 


0.9986 


0.2 


0.2 


0.0056 


0.9972 


0.4 


0.32 


0.0089 


0.9956 


0.64 


0.5 


0.0139 


0.9931 


1 


1 


0.0278 


0.9863 


2 


1.6 


0.0444 


0.9783 


3.2 


2 


0.0555 


0.9730 


4 


3,2 


0.0888 


0.9575 


6.4 


5 


0.1386 


0.9351 


10 


7.5 


0.2075 


0.9058 


15 


10 


0.2760 


0.8784 


20 



The choice of ^= 0.707 is due to the dependence of the gain on the error detector v*ihich is a 
fonction of ifae AGC and the signal to noise ratio. 

"Wlih reference to Figure 14» the next step afier timing tracking is phase tracking (step 
10 156). A block <tiagram illustrating the phase tracking portion of the variable rate modem is 
showu in Figure 17. The phase tracking comprises three main components including a phase 
detector 254, loop filter 258 and a fiequency jump 262 if Ho^ frequency is larger than a 
threshold 

The phase tracker mainly operates at the symbol rate. For higher data rates, real time 
15 constraints may pose a problem. For higjh symbol rates it may be necessary to operate the 
phase tracking at ra:tB$ slower than the symbol rate. 

The filtered I and Q data is input to amix^ 252. The output of the mixer 254 is input 
to the idiase detector 254. The QPSK phase detector performs the following: 

^„==*a^sgn(/J-/„-sgn(ej 

20 \^4iere I and Q denote in-phase and quadrature components, respectively, at a rate of once per 
symbol. 

The output of Ihe phase detector is input to a loop filter 258. The filter is presented in 
state space representation. The filter is defitied by the parameters ao and ai. The Doppler 
estimate, Dn^ is preferably also input to the filter fiom an external loop. The equations are 
25 shown below wh^-ein F„ and J*» denote the fi:equency and phase, respectively. Note that the 
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phase Pft is the same phase discussed eaxlicr denoted $ or ©, The filtesr preferably comprises 
an iijfimte impulse response (IIR) filter. The loop filter perfonns the following: 

where gn is the phase error fiximthe phase detector and ai « 1 - bj, 

If the ftequOToy Fn in the filter is large, a jump in the direct digital synthesizer (DDS) 
or other similar frequency sonrce is performed having a size tif. A typical value of JHf 
cones|K)nds to 10 Hz, 

The absolute value of flxe output of the loop filter is compared agaitist a threshold 
via a compare block 260 tihat performs the following comparison: 




If tiie expression is true then jump 4^in the opposite direction via ftequency jump 262. Note 
that the Fn and threshold values are in units of Radians/operating period. 

The data is written to the direct digital synthesizer (DDS) using frequency jump unit 
262 in accordance with tiie following 

synthesizer _ freq = synthesizer _ fhq + &f 
Note that the quantity 4f is in units of.Radians/operadng period. 

The following variables apply; 

P^^P^^^.^^F^) 
K^F„-^^sm(FJ 

It is assumed that the Doppler variable is zero. For a non-zero Doppler value one skilled in 
the art can modiigr the above equadotis to include the Doppler variable. 

Note also that the Doppler estimates are supplied from an external loop uncoupled to 
the IIR filter. 

The output of the frequency jump 262 aud the compare 260 are input to the block 
264. The output ofblock 264 fomis the second input ofthenuxer 252^ 

The transfer fimction of the IIR used in the loop filter 258 will cow be described in 
more detail, by taking the 2 transform of the two IIR equations {Fn+i, Pn-^x) and substituting 

where ^ is the phase of the input signal (linear approximation), we obtain: 
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1 - (1 + a^)2^ + (a^i + a,)^-' 
Note fhat it is coovenieot to use normalized form for a second order syst^ l.e., to defuie 

^Systran by ^si or ^ and tt^c^ i$ deteonmed by the required lock raiige, Doppler 

offsEct or output signal to noist ratio (SNR). A damping coefficient of «f * l is reoommended 
5 for reducing cycle slips as describe in Ascheid and Meyr, "Cycle Clips in PLLs: A Tutorial 
Surv^," IEEE Transactions in Communications, Vol 30, October 1982, pp. 2228 to 2241, 
For acquisition it is preferable to use a smaller damping coefBcient in the range 
J2 

^ = 0.5 to— 'Whicli yields a fester response with an overshoot of 18% for ^ = 0.5 . 
2 

After phase tracking is accomplished^ the I and Q soft decisions are output to the 
10 Viterbi decoder (not shown) (sit&p 158). Whicther signal lock has been achieved is then 
deiffittnined by comparing the BER to a predetermined teeshold (step 160). If lock has been 
achieved^ the main tracking loop is entered. If it is determined that lock has not been 
achieved than the system goes back to the acquisition (step 161)» 

A flow diagjram illustratiDg the tracking processing portion of the variable rate modem 
IS is shown in Figure IS, The various fbnctions peifimned in tracking are substantially tiie 
same as tiiose performed during the pre*tcacking stage except for the phase tracking which 
has a smaller loop bandwidth* 

The tracking stage comprises the steps of performing a matched filter (step 220), 
AGC tracking (step 222), timing tracking (stq) 224), phase tracking (step 226), I and Q soft 
20 decisions (step 22S) and determming whether lock has been achieved (st^ 230). If it is 
determined Ibat lock is achieved, the signal detection is declared (step 232). 

In operation, the main loop advances fiom stage to stage until tracking is achieved. 
At fhat point, the BER is measured by the Vxterbi decoder. If the BER indicates failure, then 
the system returns bade to acquisition. 
25 If the BER indicates lock, then the system remains in tracking until the BER indicates 

that the system is not in lock due to a large number of errors. Other indications may be used 
such as the behavior of pha:se locking loops. 

While the invention has been described with respect to a limited number of 
embodiments, it will be appreciated that' many variations, modifications and other 
30 applications of the invention may be made. 
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What is clalnz^ is: 

1 . A satellite modtm, comprisisig 
a housing; 

a antexma for receiving and trmismitdng radio fiequency (RF) signals; 

an upconverter/do^vnoonverter coupled to said antenna and adapted to upconvert RF 
signals fiom a first fiequency band to a second frequmcy band, and adapted to 
downconvert RF signals firom said second ficqnency band to said first 
frequency band; 

an RF transceiver coupled to said upconverter/downoonverter and adapted to receive 
an output transmit signal and to generate an output RF signal therefirotn, said 
RF tiaixsceiver odaptdd to receive an input RF signal £com said 
upconverter/dowaconverter and generate an ii^nit receive signal therefrom; 

an intermediate frequoacy (IP) module adapted to receive said ii^ut receive signal 
and gencgrate I and Q signals in response thereto; 

a baseband module ada^sted to lecdive said I and Q signals and to generate receive 
data ta accordance therewith; and 

a baseband xnocMe ade^ted to generate said omput transmit signal in accordance with 
a transmit data signal hspat thereto. 

2. A receiver baseband apparatus, conqprisingi 
input means adapted to receive an I and Q signal; 

a iGrst matched filt^ ade^yted to receive said I signal and generate an I filtesred output 
therefrom; 

a second matched filter adapted to receive said Q signal and gen^ate a Q filtord 

ou^ut therefrom; 
a processor programmed to perform the steps of: 

detecting the presence of signal activity as issput to said receiver baseband 
apparatus; 

acquidng said signal once it is detected; 

pre-tracking said signal once it is detected; 

tracking said signal once it is detected; 
a decoder adapted to recdive said I output signal and said Q output signal from said 

proo^sor and to generate a decoded output therefrom; 
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a deinterleav^ adapted to generate a demterleaved output in accordance vn&i said 

decoded cii^t signal ixxput thereto; 
a forward exror correction decoder adapted to gexietate output receive data in 

accordance >vitb said deixsterleaved output signal input thereto; and 
5 a controller adapted to manage and control said input means, first matched filter, 

second matched filter, said processor, said decoder^ said deinteirleaver and said 

forward error correction decoder. 

3. The sqpparatus according to claim 2, wherein said step of detecting the presence of 
signal activity comprises the steps of: 

10 performing initial automatic gain control (AGC) on said signal; 

performing signal decimation; and 
perfoiming signal detection and fiequcncy acquisition. 

4. The apparatus according to claim 2, wherein said step of acquiring said signal 
comprises the steps of: 

15 inputting said signal to a first nmtchedi&lter, 

performing a first automatic gain control (AGC) acquisition; 

perfbmiing tuning acquisition; 

inputting said signal to a second matched fQter; 

peiforming a second automatic gain control (AQC) acquisitioix; 
20 perfbnning fine frequency estimation; and 

performing phase acquisition. 

5. The apparatus according to claim 2, wliatdn said step of acquiring said signal 
comprises Hxe step of peilbrming coarse phase acquisition on said signal, wherein said step of 
performktg coarse phase acquisition oom^es the steps of: 

25 rotating vectors 2^ representing said signal into a single quadrant by an angle d^; 

wiping off said ZaC&k) modulation; 
summing the wiped off vectora Za(®k)i 

detemiimng ^ ener^ contained within a plurality of hypotheses; and 
selecting a single hypothesis firom said plurality of hypotheses having the Tpayimnm 
30 energy. 
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6. The apparatus according to claim 2, wherein said stsfp of pre-tiacking comprises the 
steps of: 

ii^nxttmg $aid signal to artiatehed filter; 
performing automatic gain control (AQC) tracking; 
per&rming timing tracking; 
pezfomiing phase tracking; 
generating I and Q soft <kcision$; and 
determining whether signal lock has been achieved. 

7. The apparatus according to claim 2, wherein said step of pre-traoking comprises the 
steps of perfozmiiig timing acquisition on K groups, each made iq) of N DFX estimates, each 
estimate calciilated from blocks of 16 symbols, said step of p^ozming timing acquisition 
compzising the steps of. 

calculating a timing estimate ti based on a DFT for 16 contigoons symbols, i = lv.>H 

theteby obtaining N DFT estimates each based on a block of symbols; 
generating a histogram of saidN DFT estimates tj; 

classifying a timing range said gtoi^ the N estimates are in based on said histogram; 
uxLvmpping said N DFT estimates and calculating their average Tu and 
unwrapping K average estimates Ti and performing a least square fit of said K 
averages so as to generate a jSnal estimate. 

8. The apparatus according to claim 1, wh^:dn said decoder comi^ses a Viterbi 
decoder. 

9. The ^aratus according to claim 1, whcrehi said upconvertcr/downconvert^ is 
adapted to upconvert an L band signal to a C or Ku band signal. 

10. The apparatus according to claim U wherein said upoonverter/downoonverter is 
adapted to downconvcit a C or Ku band signal to an L band signal 

11. A receiver baseband sq^rparatus^ comprisiog: 
input means adapted to receive an I and Q signal; 

an I matched filter adapted to receive said I signal and generate an I filtered output 
therefrom; 

a Q matched filter adapted to receive said Q signal and generate a Q filtered output 
fiierefrom; 
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a processor prograimDecl to: 

perfozm autoiaaadc gain contcol (AGC) and generate an AG control signal 
thei^fiom; 

perform timing detection and generate an A/D clock control signal tberdSrom; 

perfbnzi phase detection and generate a voltage controlled osciUator (VCO) 
control signal tberefiom; 
a decoder ad^vted to receive said I output signal and said Q output signal fiom said 

processor and to graerate a decoded output therefiom; 
a deinterleaver ad^yted to generate a deinterleaved output in accordance vdih said 

decoded output signal input dieieto; 
a forward error correcdon decoder adapted to generate output receive data in 

accordance with said deinterleaved ou^sut signal input thereto; and 
a conUoUer adapted to manage and control said input means, I matched filter, Q 

matched filter, said processor, said decoder, said deinterleaver and said 

forward error correction decoder. 

12. The apparatus according to claim 11» v^eiein said step of detecting tibe presence of 
signal activity coiapxiaes the st^s of; 

performing initial automatic gain control (AOC) on said signal; 

performing signal decimation; and 

petfimning signal detection and frequency acquisition, 

13. The £Q)paratus according to claim 11, ^vherein said step of acquiring said signal 
comprises the steps of: 

inpxitting said signal to a first matched filter, 

performing a first automatic gain control (AGC) acquisition; 

performing timing acquisitxon; 

inputting said signal to a second matched filter; 

performing a second automatic gain control (AGC) acquisition; 

performing fine fi:cqucncy estimafion; and 

performing phase acquisition. 

14. The apparatus according to claim 11, wherein said step of acquiring said signal 
comprises the step of perfoiming coarse phase acquisition on said signal, wherein said step of 
performing coarse phase acquisition cOQ^;>rises the steps of: 
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rotating vecfsors % lepresetiftiiig said signal into a singje quadrant by an angle 
iviping off said !^(ehc) moduktion; 
smnmiog the wiped off vectors zn{&jdl 

determinii^ the energy contained ^thin a phiralify of hypotheses; and 
selectix^ a smgle hypothesis &om said plurality of hypotheses having the maximum 
eneigy. 

15. The apparatas according to claim 11, wherein said step of pre-tracking comprises ^ 
steps of: 

inputting said signal to a matched filter; 
petfonning aut<nnatic gain contiol (AOC) tracking; 
performing timing tracking; 
performing phase tracking; 
generatine I and Q soft decisions; and 
determining whether ^gnal lock has been achieved. 

16. The apparatus according to claim 1 1, wherein said step of pre-tracking con^srises the 
$tq>5 of performing imaxig acquisition on K groups, each made vp of N DFT estimates, each 
estimate calculated fiom blocks of 16 symbols, said step of performing timing acquisition 
comprising die steps of: 

calculating a timing estimate t^ based on a DFT for 16 contiguous symbols, i » 1 J^, 

thereby obtaining N DFT estimstes each based on a block of symbols; 
generating a histogram of saidN DFT estimates t^; 

classifying a timing range said grot^ the N estimates are in based on said histogram; 
unwrapping saidN DFT estimates and calculating tiieu: average Tj; and 
unwrappkig K average estimates Ti and performing a least square fit of said K 
averages so as to generate a final estimate. 
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ABSTRACT 

A satellite modem M/hose baseband modules caa be implemented mostly in software 
running on a suitable processor such as a high speed digital signal processor (DSP)< The 
modem includes a RF portion and a baseband portion, The digitized output of a matched 
filter is input to the baseband processing portion that consists of software esceouting on a 
DSP. The modem comprises an antenna coupled to an upconverter/doivnconverter. The 
upconverter/doivnconverter is coupled to a RF transceiver which fimctions to receive a 
transmit signal from a baseband transmitter module and generate a signal for input to the 
upconverter/downconverter* The RF transceiver also functions to receive a downconverted 
signal and output an IF signal for input to an IF module. The IF functions to generate I and Q 
ou;lputs tibat are ii^put to a digital baseband module. The baseband module can be 
implemented as software executing on a DSP. Also included in a technique of modulation 
wipe off ^wherein the efifeets of modulation are removed. A signal detection technique 
perfomis signal detection and firequency acquisition in the presence of multiple signals ia 
accordance with a method described herembelow. A timing acquisition technique includes 
an unwrapping technique that functions to coxxect timing etiors before av^ging is 
performed. Liz^ear changes to the time estimates are performed using LSR technic^ 
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